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Abstract 

There might be a light scalar held during inflation which is not responsible for the accelerating 
inflationary expansion. Then, its quantum fluctuation is stretched during inflation. This scalar field 
could be a curvaton, if it decays at a late time. In addition, if the inflaton decay rate depends on the 
light scalar field expectation value by interactions between them, density perturbations could be 
generated by the quantum fluctuation of the light field when the inflaton decays. This is modulated 
reheating mechanism. We study curvature perturbation in models where a light scalar field does 
not only play a role of curvaton but also induce modulated reheating at the inflaton decay. We 
calculate the non-linearity parameters as well as the scalar spectral index and the tensor-to-scalar 
ratio. We find that there is a parameter region where non-linearity parameters are also significantly 
enhanced by the cancellation between the modulated effect and the curvaton contribution. For the 
simple quadratic potential model of both inflaton and curvaton, both tensor-to-scalar ratio and 
nonlinearity parameters could be simultaneously large. 

PACS numbers: 95.85.Bh, 98.80.Es, 98.80.Cq 
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I. INTRODUCTION 



Cosmic inflation solves various problems in the standard Big Bang cosmology []J and 
simultaneously provides the seed of large scale structure in our Universe from the quantum 
fluctuation of a light scalar field, e.g., inflaton field 

A single field inflation model predicts the density perturbation which is nearly scale- 
invariant and almost Gaussian. In other words, the scalar spectral index n s is close to unity 
and a non-linearity parameter /nl is much less than unity This is consistent with the 
current limit on the local type non-linearity parameter /nl from the Wilkinson Microwave 
Anisotropy Probe (WMAP) seven-year data, —10 < /nl < 74 at the 95% confidence level {4]. 

However, besides canonical single field slow-roll inflation models, there are many possible 
mechanisms to generate density perturbation. By WMAP data, the scalar spectral index 
n s has been measured with a good accuracy, while the non-linearity parameters have been 
just weakly constrained as above. The sensitivity of the Planck satellite |s| to measure non- 
Gaussianity is as good as to probe /nl of 0(1). The non-Gaussianity could be an important 
observable to discriminate between various mechanisms of density perturbation generation. 

For example, multi-field inflation models can show large non-Gaussianity with special 



conditions during inflation 
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in the radiation dominated era 1X8 



16j . preheating [l7|, or deep 



at the end of inflation 

I. The last case includes the "curvaton" scenario [13-1231 . 
A light scalar field, curvaton, has too little potential energy to drive inflationary expansion 
during inflation. At a later time when a curvaton decays, the isocurvature perturbation of the 
curvaton field becomes adiabatic or mixed with that from the inflaton field. If the curvaton 
energy density is subdominant at its decay time, the large non-Gaussianity is generated in 



general 24j. Since an inflaton also generates density perturbations, inflaton and curvaton 
contributions to density perturbation could be comparable. This mixed inflaton-curvaton 
scenario has been also studied intensively {25 ]. 

The quantum fluctuation of a subdominant light scalar field during inflation can mod- 

nn 

ulate the efficiency of reheating by the inflaton decay [26|, I27j. This makes the reheating 
a spatially inhomogeneous process. The quasi-scale invariant perturbations of this field, 
which are isocurvature modes during inflation, may be converted into the primordial cur- 
vature perturbation during this process. Large non-Gaussianity also can be induced from 



modulated reheating 28l-l3l|. For a review of modulated reheating after inflation, see for 
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example Ref . 32| . Again, in general, the density fluctuation can be generated by both infla- 
ton and modulated reheating. Mixed inflaton-modulated reheating scenario also has been 



investigated 



3l|. 



Here one may realize that a curvaton is a light scalar field and thus naturally can play 
the role of a scalar field which modulates the reheating by the inflaton decay. The inflaton 
field may have a coupling with a curvaton if that is small enough not to disturb dynamics 
of both an inflaton and a curvaton. Nevertheless this interaction modifies the decay rate 
of the inflaton field. Since the decay rate of inflaton becomes a function of the local value 
of a curvaton field cr(x), this gives rise to a perturbation in the decay rate of the inflaton 
field and thus in the reheating temperature which is responsible for the density perturbation 
after reheating. So far few attention has been paid to couplings between the inflaton and the 



curvaton 



35] . In this study, we incorporate 



compared to self-interaction of curvatons 
the modulated reheating effects in the curvaton scenario, taking the perturbation generated 
from the inflaton field also into account. 

The paper is organized as follows. After describing the model and its dynamics in sec- 
tion Ull we consider the curvature perturbation and calculate the power spectrum, the tensor- 
to-scalar ratio and non-linearity parameters in section IIHI In section [IV] we work on simple 
models and show how the modulated reheating effect by curvaton affects the density per- 
turbation and parameter space of models can be constrained. We summarize our results in 
section [V] 



II. DYNAMICS 



Inflation is driven by the potential energy of the inflaton field, <\>. We assume that 
during inflation interaction terms of inflaton with other fields are negligible. However after 
inflation, the inflaton starts oscillating around the minimum and finally, via the interaction 
terms, decays into the standard model (SM) particles, which makes the hot thermal plasma 
in the standard Big Bang cosmology. 

One of the relevant terms for the decay of inflaton field including curvaton field a could 
be given by 

C mt = A|$| V, (1) 
where $ is another light scalar field such as the SM Higgs field. Then, in the classical 
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background of curvaton field, Eq. ([I]) induces the decay of inflaton into two Higgs scalars, 
with the curvaton expectation value dependent (CD) decay width 

rTV) = ;r^V, (2) 

with being the inflaton mass at the minimum. Together with the other curvaton inde- 
pendent interactions such as 

C = O su -±-, (3) 

where Osm denotes a SM operator, induce the curvaton independent (CI) decay width of 
the inflaton, T^ 1 . The total decay rate of the inflaton is given by 

W = rJ' + rfV). (4) 

For the very light curvaton field, > with m a being curvaton mass, the curvaton starts 
to oscillate in the radiation-dominated epoch well after the inflaton decays, when the Hubble 
parameter becomes as small as m a . After the reheating by the inflaton </> is completed, the 
energy density of the radiation from inflaton decay decreases as 

p r = 3M|W^) 4 , (5) 



V a 

where a^^ is a scale factor when inflaton decays, i.e. H = T^. The energy density of the 
curvaton after the onset of the oscillations decreases as 



1 /n \ 3 

-L 2 2' Uirr1 ~ » 



with a* being the expectation value of curvaton during inflation and a ma is a scale factor 
when the curvaton start oscillation at H = m a . The curvaton decays at a later time and we 
call its decay rate r a . Whether the Universe is curvaton dominated or radiation dominated 
at the moment of curvaton decay depends on the size of T a and a*. 



III. PRIMORDIAL CURVATURE PERTURBATION 

We consider that inflaton and curvaton a fields are relevant to the density perturbation 
in the early Universe. Their vacuum fluctuations are promoted to a classical perturbation 
around the time of horizon exit. During inflation the field trajectory is dominated by the 
inflaton field and thus the inflaton perturbation becomes adiabatic and that of the curvaton 
contributes to the isocurvature mode. 
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A. Modulated reheating from an interaction with a curvaton 



Reheating of the Universe is attained from the decay of the inflaton field 0. Since the 
inflaton decay is modulated by the curvaton field a, the curvature perturbation of the 
radiation produced from the decay of inflaton has two origins. One comes from the inflaton 
field itself in the standard picture of the generation of fluctuations. The other comes from the 
light scalar field (curvaton) a during the reheating process due to the interaction between 
the inflaton and the curvaton field. Thus, as in the inflaton-modulated reheating mixed 
scenario, it is written 



1 V„ . 1 / VVsAr.o 1 / V M VVsss . vv (i 
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^ Mp V * 2M P \ VI J V * QM P \ V A 

+ Q a 5(T* + ^Q a Ml + ^QaaaSal + ■ ■ ■ , (7) 

where Q is a function of T^^/Hc calculated at a time t c which is after several oscillations 
of the inflaton but well before the time of decay of inflaton. A quantity with * is evaluated 
when the corresponding scale crosses the Hubble horizon during inflation. 

During inflation the energy density of a field is negligible and the inflation is driven by 
the inflaton field <fi alone. The slow-roll parameters during inflation are defined by 

,-^(^) 2 , ri = M^, e = M p Y^. (8) 
Using Eq. (JHJ), the curvature perturbation of radiation ( r is expressed as 

Cr = Crl + kr2 + Ur3 + " ' ' , (9) 
2 



with 



^ = u 1 /o— S< t>* + Q° Scr *' ( 10 ) 
Mpy/ 2e* 

Cr2 = ^ U ~ |r) 6fi + Q„Ml (11) 

C r3 = — ^= (-V-7T + -) + QaaaSal (12) 



e 



Mj-v 7 ^ V 2e 



B. After curvaton decay 

After inflaton decay, the energy density of radiation decreases, however the curvaton 
energy density stays the same for a while and starts to decrease when the mass of curvaton 



becomes larger than the Hubble expansion. Since the energy density of oscillating curvaton 
decreases slower than that of radiation, the curvaton becomes important well after the 
decay of inflaton. When the decay rate of curvaton T a becomes comparable to the Hubble 
expansion rate, the curvaton a decays quickly to radiation. 

After the decay of the curvaton, the remnant radiation is a mixture from the inflaton 
and curvaton decay products with different density perturbations. In this inflaton-curvaton 
mixed scenario, the curvature perturbation after the curvaton decay can be expressed as 
analytically with instant decay approximation by 
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38] 



C = Ci + 



6 



(13) 



where 



Cl=(l-i*)Crl+i2Crl, 

C 2 =(l - R)Cr2 + R(*2 + R(l -R)(S + R) (Crl - Ci) 2 , 

C 3 =(1 - R)Cr3 + R(«3 + 3i?(l -R){3 + R) (Crl " Crl) (Cr2 " Cr 2 ) 

+R(1 -R)(S + R)(-3 + 4R + 3R 2 ) (( rl - Ci) 3 , 



(14) 



with 



R 



3p 



(15) 

Here C- is given in Eq. ()9]) and ( a is the curvature perturbation from the curvaton field. R 
parametrizes the dominance of the curvaton energy density when it decays. 



For the quadratic potential of curvaton, Cr is given by : 

Cr 



1 1 

Crl + 7tCt2 + 7C<t3 H 

2 6 

2 5<t* 1 / 5(7* \ 2 2 / 5<j* x 



3 u * 



a,, 



9 V cr* 



+ 



(16) 
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Then from Eq. (fT3l) with Eqs. (1T21 - (|T6|) we obtain 



Ci 



1 - R)Q a + 

OCT* 



C-(l-«)^(l-|; 



+R(1-R)(3 + R) 

i-R ( e v 

V-7T + — 

2e e 



50* 2 8a, 

— — j== + Q a bo* - - — 
MpV2e* 3 cr* 



2i?/5a* 
3 \ er* 

2 



+3i2(l -i2)(3 + i2) 



;i-i?)g_K 3 + ^f— 

3 \ ex* 



50* 2 5a* 

— — = + QM* - - — 

MpV2e* 3 cr* 



Mj, 



+R(1 -R)(3 + R)(-3 + 4R + 3R 2 



3 \ cr* 



50* , ~ , 2 5<r* 
MpV2e* 3 cr* 



(17) 



(18) 



(19) 



C. The power spectrum 



The power spectrum of the curvature perturbation is defined by 

27T 2 

(C(*i)C(*a)> = (2vr) 3 5(k 1 + k 2 )— P c (^), 

and the perturbations of the fields at the horizon exit satisfy 

2tt 2 
k 3 
2tt 2 

k 3 

.(kJSa^fo)) = 0, 

with 



*(fci)50*(£; 2 )> = (27r) 3 5(k! + k 2 )^-^. fa} 
(M^M**)) = (2 7 r) 3 5(k 1 + k 2 )^-7V (A*} 



2^ 



(20) 



(21) 



(22) 



which is determined at around horizon exit. 

Using Eq. (1TT1) and Eqs. (1201 - (122"]) , the power spectrum of the curvature perturbation 
is given by 



with 



9 M 2 p 

p£ * [3Q^(1-R) + 2R] 2 . (24) 



9^(1 - R) 2 

Here < R < 1 defined in Eq. ( JT5l) parametrizes the contribution of the curvaton. For 
R = 1, which means that the curvaton dominates the background energy density when it 
decays, the power spectrum is given by 

n,o=(^V. (25) 

Therefore the contribution to the power spectrum from the curvaton a dynamics scales 
R 2 , while those from both the inflaton and the modulated reheating effect through the 
a field by inflaton <fi decay scale (1 — R) 2 . As one can see, in the limit of R — > 0, the 
usual curvaton contribution disappears and it corresponds to the inflaton-modulated mixed 
scenario. The opposite limit with R — > 1 corresponds to the pure curvaton scenario. Whereas 
the R parametrizes the a contribution as the curvaton compared with the other two, the a 
contribution via the modulated inflaton decay is parametrized by Q a . 

The parameter f compares the contribution to the Power spectrum from the a field to that 
from inflaton 0. In the limit of f — > 0, the Power spectrum comes solely from the inflaton, 
while in the limit of f 3> 1 the a field contributes dominantly through the modulated effect 
and/or the curvaton effect with the inflaton effect suppressed. 

One should notice that because the latter two contributions come from the single same 
source a, there is a cross term of two. This cancellation between the modulation effect and 
the curvaton effect makes non-trivial features in the Power spectrum. Both contributions 
may cancel each other when 3Q a (j* ~ —2R/ (1 — R) and the inflaton contribution dominates. 
To quantify the amount of the cancellation, we define 5 as 

6 = 1+ 3 ^( 1 ~- R ) (26) 
2R 

This 6 is a measure of fine tuning of the cancellation and becomes 5 = for the exact 
cancellation. The contour plot of 6 is shown in figure [TJ With this, f can be written as 

= - 8M P e * R2 6 \ (27) 
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From Eq. (1231) . ignoring the negligible contribution from the curvature of the curveton 
potential V aa , the scalar spectral index n s is given by 

^- ls gL *. + Z±±*k, (28 ) 
dmk 1 + r 



0.0 0.2 0.4 0.6 0.8 1.0 
R 



FIG. 1: The contour plot of 5 defined in Eq. (|26p . The cancellation happens along the line of 
5 = 0. 

which we normalize to be 0.97 through in our analysis. The tensor-to-scalar ratio is given 
by 

16e 

TT 55 n = d-Rm + ry (29) 

where we used Vt = 8(H*/2ir) 2 . As you can see here, for small f, the observational limit 
r T < 0.36 constrains the value of R to be R < 1 — a/166^/0.36 ~ 0.53. 

In figure [2j we show the contour plot of f for cr* = 0.05 Mp (which we will later call 
Case B) with e* ~ 0.005 and 77* = 0. There is a cancellation between curvaton effects and 
modulated effects around the dashed line (blue) which connect (R,Q a a*) = (0.6,-1) and 
(0, 0), where f vanishes. In this small f limit the inflaton contribution dominates the power 
spectrum. In the opposite region with a large f, the a field dominates the Power spectrum. 
For different values of cr*, the magnitude scales as <r~ 2 , since e* does not change much. 

In the upper left panel of figures |3] - El we show the power spectrum in the plane of 
(R,Q a a*) plane in the upper-left window for each cases with different cr* ; 10~ 3 Mp (Case 
A), 0.05 Mp (Case B), 0.5 Mp (Case C). We have normalized the amplitude of the power 
spectrum by the value at the pure curvaton limit of R — 1, namely V^/V^q. The slow-roll 
parameter e* is calculated from Eq. ( j28j) considering n s = 0.97 and we used for simplicity 




FIG. 2: [Left window] : The contour plot of f for Case B, <7* = 0.05Mp. For the other cases 
the magnitudes are scaled by o\. Along the blue dashed line f = 5 = . [Right window] : 
The tensor-to-scalar ratio for Case B. The blue shaded region is ruled out using the constraint 
r-r < 0.36. 5 = along the blue dashed line. 

?7* = = 0. We assume that the correct observational value can be attained using the 
residual parameter of H*. 

The three cases in figures [3] - [5] show the following features: 

• Case A : cr* = 10~ 3 Mp and f ^> 1 in most of the region. 

Modulated reheating and curvaton contributions are dominant and inflaton contribu- 
tion to the curvature perturbation is subdominant, 

• Case B : cr* = 0.05 Mp and f ~ 1 as shown in figure [2j 
All three contributions are effective, 

• Case C : a* = 0.5 Mp and f 1 in most of the region. 

Inflaton contribution is dominant and it scales as (1 — R) 2 . R < 0.45 is allowed from 
the constraint on the tensor-to-scalar ratio. 

In the Case A (upper- left window in figure [3]), 5a contribution is dominant in the overall 
region and the inflaton contribution is subdominant. In the cancellation region, the relative 

10 



p d p w /nl 




R R 

FIG. 3: The contour plots on R — Q a a^ plane, of the power spectrum normalized by the 
value at R = 1 (upper- left), /nl (upper-right), tnl (lower- left) and #nl (lower-right) for Case 
A, cr* = 1CT 3 M P . We put r?* = £* = 0. The red shaded region corresponds to too large non- 
Gaussianity to be consistent with the observation, —10 < /nl < 73, tnl < 10 4 , and |<7nl| < 10 5 . 
Along the blue dashed line f = the cancellation happens. 

contribution from 5a decreases. Of course, by increasing H*, the desired amplitude of 
can be recovered. The tensor-to-scalar ratio is always much smaller than 0.01 in the shown 
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FIG. 4: The contour plots on R — Q a a^ plane, of the power spectrum normalized by the 
value at R = 1 (upper- left), /nl (upper-right), tnl (lower-left) and ^nl (lower-right) for Case B, 
a ^ = 0.05 Mp. We put 77* = = 0. The blue shaded region corresponds to too large tensor-to- 
scalar ratio to be consistent with the observation, rj- < 0.36. Along the blue dashed line 5 = 0. 

parameter range. 

In the Case B (upper- left window in figure H]), the inflaton contribution is effective for 
cancellation region and modulation and/or curvaton effects are effective in the rest. The 
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FIG. 5: The contour plots on R — Q a a^ plane, of the power spectrum normalized by the 
value at R = 1 (upper-left), /nl (upper-right), tnl (lower- left) and (?nl (lower-right) for Case C, 
o"* = 0.5 Mp. We put ?7* = <!;* = 0. The blue shaded region corresponds to too large tensor-to-scalar 
ratio to be consistent with the observation, < 0.36. Along the blue dashed line 5 = 0. 

cancellation between the modulated reheating and the curvaton appears around R ~ 0.7 for 
negative Q a <J*. The tensor-to-scalar ratio is comparable to the observation constraint and 
a region around (R, Qo-a*) ~ (0.6, —1) is ruled out as shown in the right window of figure EJ 
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In the Case C (upper-left window in figure [5]), the density perturbation dominantly comes 
from the inflaton field but the magnitude changes due to the effect of curvaton. For large 
R the magnitude decreases. A portion of parameter space is excluded because of too large 
tensor-to-scalar ratio. 



D. Non-Gaussianity 

When the curvaton field dominates the energy density when it decay, i.e. R ~ 1, the 
curvature perturbation is dominated by the pure curvaton and the non-Gaussianity is sup- 
pressed. However in the other case < R <C 1, there is a possibility to get large non- 
Gaussianity. 

The bispectrum is given by 

(CfeCfeCs,) = (^) 3 B c (h, k 2 , k 3 )5(h + k 2 + k 3 ), (30) 
and the dimensionless non-linearity parameter for the bispectrum, /nl, is defined by 

B c (h, k 2 , k 3 ) = ^hL[Pc(h)P c (k 2 ) + P c {k 2 )P c {k 3 ) + P c (fc 3 )P f (A*)]. (31) 
From the trispectrum Tq given by 

(C^C^c = (2tt) 3 T c (A ;i , k 2 , k 3 , k 4 )5(h + h + k 3 + h), (32) 



the dimensionless non-linearity parameters tnl and #nl are defined as [39] 



T c (k u k 2 , k 3 , k 4 ) ^^[P^k^P^k^P^h) + 11 perms] 
54 

+ —g^P^P^P^h) + 3 perms]. 



(33) 
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, .4 if . For instance, Smidt et. 
.3 [40|. 



The current bounds have been derived by several groups 
al. reported as -7.4 < # N l10 -5 < 8.2 and -0.6 < r NL 10~ 4 < 3. 

In figure [3], we show the contours of the non-linearity parameters /nl (upper-right win- 
dow), Tnl (lower- left window) and #nl (lower-right window) for the Case A. For the contours 
we have taken account the relations, Q aa al = Q a <T*(6Q a <T* + l) and Q aa o&l = ^Qo&*Qoa&li 
which are motivated from the specific example we will show in the next section. We can see 
clearly the enhancement of non-linearity parameters in the cancellation region of modulated 
reheating and curvaton. The pure curvaton limit is recovered along the line of Q a a* = 0. 
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In figure H] for the Case B, we can see the enhancement of the non-linearity parameters 
along the cancellation region. The difference from the Case A is that now the inflaton 
becomes more important and diminishes the non-Gaussinaity. 

In both cases of A and B, the large non-Gaussianity is dominantly due to 5a, while the 
power spectrum comes from both depending on f , as explained in the Appendix. The non- 
linearity parameters can be enhanced around the cancellation region (5 < 1 and 5<f). In 
this region, a large /nl comes dominantly from ( 2 ,aa as defined in Eq. flA2j) of the Appendix 
and estimated to be 

, 5 r 2 19(l-R)Q aa a 2 ;-6R + R(l-R)(3 + R){Q a a*-2) 2 
/NL " 6 (ITW 4^ • (34) 

Therefore small 5 (large cancellation due to fine-tuning) can induce larger /nl- However 
note that f is proportional to 5 2 (so r 2 oc £ 4 ), thus too small 5 makes r becomes smaller 
than 5 itself and reduces /nl- One can see this behavior clearly in the upper right figure of 
figure HJ /nl decreases when approaching the cancellation line (blue dashed line). 

In the same region tnl is also dominated by £2,0-0- term and approximately the squared of 
/nl, 

5 /nl) 2 . (35) 

Large #nl is possible in the same cancellation region dominated by C3,ow term and given by 
25 r 3 11 



9nl 



54 (1 + f) 3 S 3 8R 3 

x [9(1 - R)Q <J<ja al + 9/2(1 -R){3 + R)(3Q a a* - 2)(3Q aa a 2 
+ 36R + R(l -R)(3 + R){-3 + AR + 3R 2 )(3Q a a* - 2) 3 ]. 



(36) 



As you can see here, (?nl has different sign in the opposite side of the cancellation line due 
to the odd exponent '3' of 5. 

Here one can see that there are two points in both sides of 3Q a a* ~ — 2i?/(l — R) those 
give the same values of and /nl- This means that the measurements of only and 
/nl can not determine R and Q a a* uniquely. However, this degeneracy can be resolved by 
measuring #nl because one predicts <?nl to be positive and the other does it to be negative. 

In figure [5] for the Case C, the region of large non-linearity parameters are in the region 
of R > 0.5 but which is excluded out by large r^. 
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IV. A SIMPLE MODEL 



As we have seen in the previous sections, if a light scalar curvaton field a has an in- 
teraction with inflaton field <ft, the fluctuation of curvaton field 5a can modulate reheating 
through the decay of inflaton and can affect the curvature perturbation besides the usual 
curvaton mechanism. The resultant power spectrum and the non-linearity parameters can 
be considerably affected. In this section, we examine these effects in a simple model. 

The inflaton and the curvaton a can have the following interactions in the scalar 
potential 



as well as an interaction with another scalar field as given in Eq. (TfJ. It is also possible that 
the inflaton has interaction with other fields independent from the curvaton. We consider 
that the inflation is dominantly driven by a single inflaton field with its quadratic mass term 
potential by assuming ^> m a and 



We consider cases of vanishing \ ( j )CT in subsection IIV Al and IIVB| and mention the effect of 
nonvanishing A^ CT in subsection IIV CI Its interactions with other fields are important during 
reheating or later. Therefore during inflation, field equations are reduced to 



V[(f>, cr] = -ml4> 2 + -mla 2 + -A^0V, 



(37) 



ml > X^l- 



(38) 




(39) 



(40) 



under the slow-roll condition 




V = 



M 2 P V H 
V 



«1, 



(41) 



and the inflaton- domination condition 



1 2/2^1 22 
» 2 m a a ■ 



(42) 



By solving field equations, we obtain 



^ = (4A inf + 2)M|, 



(43) 
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FIG. 6: [Left window] : The contour plot of of the Model I in the unit of log 10 fig^l 
where V( has been fixed by the observed value V( = 2.44 x 1CP 9 . In the left side of the red line 
radiation is dominated at the time of curvaton decay and in the right side curvaton dominates. The 
cancellation between curvaton and modulated reheating occurs along the blue line where r = 0. 
[Right window] : The contour plot of r. The red and blue lines are the same as in the left 
window. 

where iV in f is the number of e-fold at horizon exit from the end of inflation. The power 
spectrum is 

p < - ww (2Ar "' + 1)2 + { (1 - R)Q '°' + t}' + 1} ' (44) 

The inflaton mass controls the amplitude of the density perturbation. For the observed 
V( we can estimate m^. 

It is natural to assume that inflaton decays at its oscillating stag e after a while, so that 



r '<f,(<j) I 'm^ is very small. In this case Q is well approximated by 30] 



Q-l^mi (45) 
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FIG. 7: The contour plot of for the Model I. The blue line shows f = 0. 



The derivatives are expressed as 

Qaacr = ~ | ^ - + | . (48) 



A. Model I : Inflaton decays through only the coupling in Eq. ([!]) 



First, let us consider the case that the inflaton decays through only the coupling in 
Eq. (CD). For this case, we obtain 

(Q*CT*, Q^l QaaaVl) = ( -§, §, -§ ) . (49) 

The parameter defined in Eq. (|T5|) at the time of curvaton decay is can be obtained using 
Eqs. © - (ED by 

R = al,Ml (50) 
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with 



a 



a r CT 




) 



1/2 



1/3 



( 



) 



2/3 



radiation — dominated 



a— dominated 



(51) 



Here r CT is the decay rate of the curvaton field and a mCT /ar CT is expressed in different ways 
depending on whether it is radiation-dominated or curvaton-dominated when the curvaton 
decays at H = T a . 

In this case, there are only three parameters m^a, and T a /m a , since Q a a* and so on 
are completely fixed. As mentioned above, from the normalization of the power spectrum, 
7710 can be expressed by the other two parameters as shown in the left window of figure [61 
The largest value to m$ in figure [6] seems to be 10~ 5,2 Mp which is the same as the results in 
the quadratic chaotic inflation. This is because the region corresponds to the cancellation, 
3Qo-o"* — —2/2/(1 — R) line (blue line and R = 1/3 in this case), where f ~ and the 
contribution from a field cancels, and the dominant contribution comes from field. On 
the other hand, a contribution is not negligible in other regions and hence a smaller is 
needed to produce the observable V$. In the right window of figure [6], we showed a contour 
plot of f. We note that the condition for a negligible \^ a , Eq. (158]) . can be rewritten as 



In figure we showed the contour plot of tt- The present bound on the tt does not 
constrain this model, however future bound can rule out the parameter range around the 
cancellation region along the blue line. 

In figure [HI the non-linearity parameters /nl (left) and #nl (right) are shown. 

B. Model II : Inflaton decays through the coupling in Eq. (fTj) and others 

Next, let us consider the case that the inflaton has a nonvanishing a independent decay 
modes and the total decay width is given by Eq. (J4j). For this case, we obtain 



with Br = /IV Obviously, the Br —> 1 limit reduces to the Model I in the previous 
subsection. In this case, we have four parameters m^,,(T„r ' a / m a and Br. Again, can 
be used for the normalization of the power spectrum. In figure [91 we show the contour plot 




(52) 
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FIG. 8: [Left window] : The contour plot of /nl m (^*/Mp,T a /m a ) plain for the Model I. The 
red shaded region corresponds to too large /nl to be consistent with the observation. [Right 
window] : The contour of #nl for the Model I. 

of /nl for Br = 0.3 (left) and 0.1 (right). As we have seen, the cancellation happens at 
3Q CT a* ~ -2/2/(1 - R) or R ~ Br/ (2 + Br). 

C. Model III : The effective inflaton mass with the coupling in Eq. ([3]) 

Finally, let us consider the case of nonvanishing A^ CT . Then, through just only the coupling 
([3]), in other words even without direct coupling (fITh the inflaton decay width T ~ M^/Mp 
with the effective inflaton mass M| = m| + A^of which is a dependent. 

For this case, we obtain 

(Qa(J*, Qv*°l, Qcca°l) = ( -\Fr, \Fr(2Fr - 1), Fr 2 (3 - 4Fr) ) , (54) 

with Fr = X^al/M^. The qualitative behavior of this model is the same as that of the 
Model II up to numerical factors, by replacing Br with Fr. 
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FIG. 9: [Left window] : The contour plot of /nl with Br = 0.3 on (a^/Mp, T a /m a ) plane for the 
Model II. The red shaded region corresponds to too large /nl to be consistent with the observation. 
[Right window] : The same as left window but with Br = 0.1. 



V. CONCLUSION 



We have studied the case where a light scalar field a induces the modulated reheating 
by the inflaton decay and also acts as the curvaton by its late time decay in the presence of 
the curvature perturbation generated from the inflaton field itself. In fact, the coupling in 
Eq. ([T]) is possible from the gauge invariance of the SM 1 , provided both and a are gauge 
singlet as is often assumed to preserve the flatness of the potential. 

When u field contributes to both inducing the modulated reheating by the inflaton decay 
and generating the density perturbation as the curvaton, there could be a cancellation 
between two contributions along the line 3<5o-cr* — —2/2/(1 — R). Around this cancellation 
region, 8a contribution to the power spectrum of the density perturbation is subdominant 
and the inflaton contribution becomes dominant. Near such a parameter region, mostly the 
middle range of R and a negative Q a cr*, non-linearity parameters tend to be large because 

1 Of course, it is also possible that those terms are forbidden by additional symmetry such as a certain 
Z2-parity. 
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of the cancellation between the modulated reheating and the curvaton originated from the 
same field a. In this sense, this cancellation is a kind of mechanisms to generate a large 
non- Gaussianity 

As specific models, we have also studied a quadratic inflation and curvaton model and 
demonstrated how the parameter space of a given inflaton and curvaton model would be 
constrained by taking the interaction between them into account. The measurement of non- 
linearity and tensor-to-scalar ratio may probe the strength of (non-)interaction between the 
inflaton and the curvaton. 
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Appendix A: Large non- Gaussianity 

The curvature perturbation £ can be decomposed of the contributions from each field 
perturbations order by order as 

C = Ci + ^C 2 + ^C3 + ..., (Ai) 

with 

C2 <2 M m 2 + 2C 2 ,^(50)(H + C 2 ,^(H 2 , , . 

(A2) 

C 3 =C3,^(50) 2 + 3C 3 ,^(^) 2 (5a) + 3(3,^(50) (H 2 + C,^{5af. 

Here we considered two-field case for simplicity and we suppressed * which denotes the 
value at horizon exit, i.e. 5<p = 5(f)* and 5a = 5a*. For the model of modulated reheating by 
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curvaton, we can read each component from Eq. (fl9|) . 
The power spectrum is given from Eq. f l20l) by 

where we used Eq. (1221) and f = Ci^/Ci,^- From the definitions in Eq. fl3Tj) and Eq. ( 133]) . the 
non-Gaussianity parameters are given by 



/] 



NL 



6 (CL + CL 



2 N2 



,ac 

TNL = , > 2 \3 7 ( A4 ) 



#NL 



(Cl,<^ + Cl,cr) 3 

25 2~^a 6 c Cl,aCl,&Cl,cC3 



54 (d 2 , + d 2 .) 3 ' 
Here a, 6, c denotes and a. 

First we will investigate the condition for large /nl in the curvaton-modulated scenario. 
Using Eq. (119j) . the denominator of /nl in Eq. ( 1A4I) becomes 



(Cl,(4 + Cl,(T 



2 \2 



1 



2M 2 e* 



(1 - R)\l + f 



;A2 



;i + f 



(A5) 



where in the last equation we dropped 0(1) coefficient, R and Q CT cx. In the same way, the 
numerators are 

tl&M =7T7Tt( 1 - R?[R{R + 3) + 2e* - 77 1 



2&XvC2,*t = QM 2 1 £ " ^) 2 ( 3 + ^)[3Q^(1 - i?) + 2i?][3Qa - 2] ~ 

CvCa.™ =^[3Q <T <7(1 - i?) + 2 J R] 2 [9Q a(7 a 2 (l - R) - 6R + R(l -R){3 + R)(3Q a a - 2f 
' 1 f 2 

(A6) 

Therefore we find that large /nl can arise from the last term Ci^.ao to give 

f 2 1 



It is easy to see that 



/nl~^ for f>l, and 8 < 1. (A8) 
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and 

/nl~tt for r < 1, and 5 < f. (A9) 

However f and 5 are not independent variables and f is proportional to 5 2 in our case as 
in Eq. (12 7p . Therefore for a given parameters, /nl vanishes in the limit of 5 — > when f 
becomes vanishing too. The large /nl can arise whenl 5 is smaller than 1 and smaller than 
r. 

Similarly, large tnl and <?nl are obtained in the same region as that of /nl when £ 2 ,oo 
and £3,0-0-0- term dominates and are estimated to be 

TNL ~(ITW and ^ NL ~(TTW (A10) 

We find that #nl is smaller than r NL in the large non-Gaussianity region. 
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